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Abstract. We have used the LONGSP spectrometer on 
the 1.5-m TIRGO telescope to obtain long slit spectra in 
the J, H, and K wavelength bands towards two positions 
along the Orion bar. These data have been supplemented 
with images made using the ARNICA camera mounted on 
TIRGO as well as with an ESO NTT observation carried 
out by Dr A. Moorwood. We detect a variety of transitions 
of hydrogen, helium, OI, Fell, Felll, and H2. From our 
molecular hydrogen data, we conclude that densities are 
moderate (3 — 6 x 10 4 cm~ 3 ) in the layer responsible for 
the molecular hydrogen emission and give no evidence for 
the presence of dense neutral clumps. We also find that 
the molecular hydrogen bar is likely to be tilted by ~10 
degrees relative to the line of sight. We discuss the relative 
merits of several models of the structure of the bar and 
conclude that it may be split into two structures separated 
by 0.2-0.3 parsec along the line of sight. It also seems likely 
to us that in both structures, density increases along a line 
perpendicular to the ionization front which penetrates into 
the neutral gas. 

We have used the 1.317/im OI line to estimate the 
FUV radiation field incident at the ionization front and 
find values of 1 — 3 x 10 4 greater than the average inter- 
stellar field. From [Fell] line measurements, we conclude 
that the electron density in the ionized layer associated 
with the ionization front is of order 10 4 cm~ 3 . Finally, our 
analysis of the helium and hydrogen recombination lines 
implies essential coincidence of the helium and hydrogen 
Stromgren spheres. 

Key words: interstellar medium: HII regions - interstel- 
lar medium : Orion - IR: interstellar medium : lines and 
bands 



1. Introduction 

The properties of the Orion nebula are the starting point 
for many of our ideas on high mass stars and their in- 
teractions with the environment. A symposium held in 



1981 (Glassgold et al. 1982) summarises much of what was 
known at that time. More recent work has been reviewed 
by Genzel & Stutzki (1989). In general, the aim has been 
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to understand the ionization structure and dynamical evo- 
lution of the nebula. In recent years, much attention has 
been paid to the hot neutral gas adjacent to the ionization 
front known as a PDR or Photon Dominated Region. 

Work on the ionized nebula has tended to concen- 
trate upon determinations of the ioniz ation structure 
and e lemen tal abundances (see Peimbert 1982 , Simpson 
et al. 1986 ). More recent optical and infr ared studie s 
have been carri ed out by Osterbrock et al. (1990 , 1992 ), 
Bal dwin et al. fll99lD , Pei mbert et al. ( |l992[ ), Po gge e t 
al. (|1992| ), DePoy & Pogge ( |1994[ ), Bautista et al. ( |l995[ ), 
Rubin et al. ( 1993| ), and Rodriguez ( 1996) . A review of the 
results has been made by Peimbert (1993) and discussions 
of the methods employed are given by Mathis dl995| ) and 
by Peimbert ( 1995 ). These studies in general indicate that 
the major fraction of elements such as C,N,0,S are in the 
gas phase within the ionized nebula whereas species such 
as Si and Fe appear to be depleted by roughly an order 
of magnitude relative to abundances either in the Sun or 
nearby B stars. 

Radio work on Orion has revealed an immense variety 
of structures in the emissions of the ionized gas ( e.g. Fclli 
et al. 1993, Yusuf-Zadeh 1990] ). Particularly striking is the 
bar-like structure situated roughly 2 arc minutes (0.25 
parsec) to the south-east of the Trapezium stars which is 
the subject of this article. "The Bar" is also observed in 
molecular line emission (see below) and clearly marks an 
ionization front where Lyman continuum photons from the 
06 star <d 1 G Ori are absorbed. The dynamical behavior 
of the ionized gas can be studied in radio recombination 
lines (Panko nin et al. |l979| , Wils on fc Jager |l987t Wil- 
son & Filges 199C| , Wilson et al. 1997 ) from which one 
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concludes that much of the ionized material in Orion is 
streaming towards the observer. 

Infared studies of HII regions sample not only the ion- 
ized gas but also the adjacent neutral material or PDR. A 
recent review of the pr opert ies of these regions is that of 
Hollcn bach fc Tielens ( 1997) (se e also the discussions of 
Genzel 1992 , and Walmsley 1997 ). Modelli ng stu dies have 
been carried out by Tielens & Hollenbach (1985), Hollcn- 
bach et al. fll99l|), Sternberg & D algarn o ( pJ89| ; |1995| ), 
Fuente e t al. (|1993"| ), Jansen et al. ( |l995a|,b), Bertoldi & 
Draine fll996[) and Draine & Bertoldi ([l996|). Much of 



this activity has centred on attempts to understand the 
properties of "The Bar" mentioned above. Recent obser- 
vational studies using a variety of molecular tracers have 
been carried out by Tielens et al. ( 1993]), b y Tauber et 
al. ( [1994[ |1995[ ) by H ogerheijde et al. ( |1995| ), and by van 
der Werf et al. (1996). These show a stratification along 



the direction of the perpendicular to the bar in the plane 
of the sky. This is in the sense expected for gas heated 
by the Trapezium stars and consistent according to the 
models with attenuation by a gas of density 5 10 4 cm~ 3 . 
However, the data also seem to show that the gas in the 
bar is far from homogeneous and that clumps of density 
as high as 10 6 cm~ 3 are embedded in the filament. Such 
high density condensations presumably either have been 
or will be soon overrun by the ionization front and will 
give rise to dense ionized globules within the HII region 
(see e.g Lizano et al. 1996, Dyson et al. 1995). Under- 



standing the characteristics of such high density clumps 
may thus be of critical importance for the evolution of the 
HII region. 

One of the most useful tracers of PDR's has turned 
out to be the near infrared lines of molecular hydrogen. 
For example, van der Werf et al. used the FAST camera 
on the ESO/MPI 2.2 m telescope to image the H2 v = 
1 -> S(l) (2.122 H and v = 2 -> 1 S(l) (2.248 fxm) 
lines towards the bar with 1.5" resolution. These show 
that the transition from atomic to molecular hydrogen 
in the bar occurs 15 arc seconds (0.03 pc at a distance 
of 450 pc) to the SE of the ionization front (i.e. away 
from the ionizing stars). Van der Werf et al. also find that 
the ratio R12 of the intensities of the 1 — > and 2 — » 
1 lines varies between a value of 8.1 ± 0.7 at the peak 
of the H2 emission to a value of 3.4 ± 1.9 30" from the 
ionization front on the side shielded from the radiation 
of the Trapezium stars. The latter value is characteristic 
of UV-pumped fluorescent emission in a low density gas 
(Sternberg and Dalgarno 1989). 

The present study had as its aim to obtain near IR 
spectra of the gas in the vicinity of the bar in order to 
verify and extend understanding of the physical conditions 
on both sides of the ionization front. We were partly moti- 
vated by the idea that there is a link between the ionized 
and PDR components in that the former is mainly sen- 
sitive to the radiation just shortward of 912 A while the 
latter is basically a measure of the radiation longward of 



this limit (see Bertoldi & Draine 1996 for a discussion). 
It is thus of considerable interest to compare the two us- 
ing the same instrument. We therefore used the TIRGO 
telescope on the Gornergrat (Switzerland) to obtain slit 
spectra in the J, H, and K bands at 3 positions in the 
vicinity of the bar, shown in Fig. [I]. As supplementary in- 
formation, we also made use of unpublished observations 
carried out using the IRSPEC spectrometer on the ESO 
NTT telescope by Dr A.Moorwood. 

We summarize in the next section the techniques used 
for observations and data reductions. The results are pre- 
sented in Sect. 3 and discussed in Sect. 4. We summarize 
our conclusions in Sect. |[ 

2. Observations 

2.1. ARNICA Observations 

The Orion Bar was observed during two observing runs in 
January 1996 and February 1997 using ARNICA (ARcetri 
Near Infrared CAmera) mounted on the 1.5m TIRGOEl 
telescope. ARNICA is equipped with a 256x256 NIC- 
MOSS array, the pixel size with the optics used at TIRGO 
is 0.96"; for a complete description of the instrument 
and of its p erformance, see Lisi et al. ( 1996 ) and Hunt 
et al. ( 1996| ). The Bar was imaged in the three J, H, 
and K broad band filters (centered at 1.25, 1.65, and 
2.2 /im, respectively) and in the Br7 narrow band filter 
(A = 2.166 /mi, AA/A - 1%, Vanzi et al. |l997b| ). The 
seeing was approximately 2-3" and the observed field was 
approximately ~ 4'. 5 x 4'. 5, covering all the Bar region. 
Data reduction was carried out using the IRAF la and AR- 



NICA a (Hunt et al. 1994) software packages. Photomet- 
ric calibration in the J, H and K bands was performed 
by observing photometric standard stars from the list of 



Hunt et al. (1997) ; the calibration accuracy is estimated 
to be ~ 5%. The Br7 image was continuum subtracted 
and calibrated using the K band image. We show in Fig. [I] 
(left panel) an image obtained combining the J, H, and K 
images and (right panel) the continuum subtracted Br7 
image. 

2.2. LONGSP Observations 

J (1.25/im), H(1.65 /im), and K (2.2/im) band spectra of 
the Orion Bar were obtained using the LonGSp (Longslit 
Gornergrat Spectrometer) spectrometer mounted at the 
Cassegrain focus on the TIRGO telescope. The spectrom- 
eter is equipped with cooled reflective optics and a grat- 



1 The TIRGO telescope is operated by the C.A.I. S.M.I.- 
C.N.R Firenze, Italy 

2 IRAF is made available to the astronomical community by 
the National Optical Astronomy Observatories, which are op- 
erated by AURA, Inc., under contract with the U.S. National 
Science Foundation 

3 A descriptio n of ARNICA can be obtained from the Arcetri 
Observatory at |itp : //150 . 217 . 20 . 1 /pub/ arnica/ 
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Fig. 1. On the left: composite J, H, and K image of the Orion Bar region showing the slit positions used for the LONGSP and 
IRSPEC observations. Coordinates are offsets in right ascension and declination relative to the position of the star 2 A Ori 
(R.A(1950)= 5 h 32 m 55. s 5 , Dec(1950)= -5° 26' 51"). We show our three slit positions (l,2,and 3) as well as the four positions 
for which we tabulate line intensities (A,B,C,CS). On the right: Br7 image of the Orion Bar region. 



ing in Littrow configuration. The detector is a 256x256 
engineering grade NICMOS3 array (for detector perfor- 
mances see Vanzi et al. 1995). The pixel sizes are 11.5 A 



(first order) and l'/73 in the dispersion and slit directions, 
respectively. LONGSP operates in the range 0.9-2.5 yum 
achieving a spectral resolution at first order of R ~ 550 
in J, 700 in H and 950 in K. For a more comprehensive 



description of the instrument, refer to Vanzi et al. (1997a) 



Observations were conducted in two runs in January 
and March 1996 under non-photometric conditions. The 
slit used had dimensions 3'.' 5x70" and was oriented N-S. 
The seeing during the observations was in the range 2"-4". 
The Orion Bar was observed at three slit positions labeled 
as 1, 2, 3 and shown in Fig. [j] superimposed on a NIR 
image obtained by co-adding the J, H and K ARNICA 
observations discussed in the previous section. Position 1 
and 2 were chosen in order to study the variation of line 
intensities along a cut encompassing all the bar. Position 
3 was subsequently chosen to be coincident with the CS 
peak discovered by van der Werf et al. at R.A.= 5 32 m 
58.5 s and Dec.= -5° 26' 25" (B1950.0). This high density 
(10 6 cm -3 ) clump appears to be illuminated directly from 
the Trapezium and we thought it useful to examine di- 
rectly the relative variations in line intensities across the 
clump. The center of the slits were offset by -35", -23" 
(Pos. 1) -35", 23"(Pos. 2) and 42", 26" (Pos. 3) in R.A. 
and Dec. with respect to the star 8 2 A Ori. 



At each grating position we performed 5 ABBA cycles 
(A=on source, B=on sky) with an on-chip integration time 
of 60 sec, for a total of 10 min integration on source. At 
the beginning or at the end of the five cycles on the object 
we performed 1 ABBA cycle on the 06 star BS 1895 (S 1 C 
Ori). 

Data reduction was performed with the ESO package 
MIDAS, within the context IRSPEC, modified to take into 
account LonGSp instrumental characteristics. The frames 
were corrected for bad pixels, flat-fielded, sky subtracted 
and wavelength calibrated using the OH sky lines present 
on all the frames (Oliva & Origlia 1992 ). After a direct 
subtraction, sky removal was optimized by minimizing the 
standard deviation in selected areas where the OH sky 
lines were poorly subtracted but no object emission was 
present. The wavelength calibration was performed to bet- 
ter than 1/5 of a pixel (~2A). 

The spectra were then corrected for telluric absorption 
by dividing by the featureless spectrum of : C Ori. For 
more details on LonGSp data reduction, see Vanzi et al. 



1997a 



Flux calibration of the spectra in the J, H, and K 
bands was achieved by rescaling the observed flux dis- 
tribution along the slit to match that obtained from the 
ARNICA images at the positions of the slits. We consider 
such calibration accurate to ~20% when comparing the 
fluxes of lines measured in two different bands. Indeed, the 
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Table 1. Spectra in the J, H and K bands 
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comparison between the flux distributions of H2(1,0)S(1) 
A2.12/im in our observations and in those by van der Werf 
et al. shows only a 10% discrepancy in the absolute flux 
level. 

2.3. IRSPEC Observations 



2.08 x 10" 3 erg cm^s" 1 sr~ x 



C ) 0.61 x 10~ 3 erg cnrV 1 sr" 
CS) 0.50 x 10~ 3 erg cm'V 1 sr" 



IRSPEC (Moorwood et al. |l99l|) observations of the bar 
using the ESO NTT telescope were carried out in 1991 by 
Dr A. Moorwood. The detector was a SBRC 62x58 InSb 
array with pixels of ~5A (H band ) along the dispersion 
and 2'.'2 along the slit direction. The slit, 4'.'4xl20" in 
size, was oriented NE-SW as shown in Fig 0. The data 
were uncalibrated and in this paper, we merely make use 
of the profiles of line intensity along the slit. 

3. Results 

In Figs. ^, ^, ^, we show sample spectra averaged over 
three portions of slit positions 1 and 2 (see Fig. p. Based 
on the profiles of line intensity along the N-S direction, 
we decided to divide the combined slit into three sections 
which we named A (a 28 pixel section to the north), B (14 
pixels in a central region), and C (28 pixels to the south). 
These sections are displayed in the left panel of Fig. 1. The 
hydrogen and helium recombination lines peak in position 
A and become weaker to the south whereas the molecular 
hydrogen lines become stronger and reach their maximum 
intensity in position C. However, there is clearly ionized 
gas towards the region C and vice- versa. In Figs. |[ |^, |[ 
we show also the spectra summed over 40 pixels in slit 
position 3 (the "CS peak"). 

In Table 1, we give intensities corresponding to the 
spectra shown in Figs. ||, H||. They have been averaged 
over the respective apertures. Line intensities, in each re- 
gion, are given relative to Br7 (put equal to 100). Typical 
uncertainties vary from ~ 10% for strong lines (/ > 10) 
to ~ 20% for lines with 1 < 7 < 10 and about 50% for 
the others. When comparing lines in two different bands, 
a 20% error due to spectrophotometric calibration must 
also be taken into account. 

The line intensities can be corrected for reddening 
using the Cardelli et al. (1989) prescription A\/Ay = 
0.48A^ 1,61 , where we have adopted i?=5.5, as appropriate 
for the Orion region. ^From the ratio Pa/3/Br7 we derive 
Ay ~ 2 mag (Sect. 3. 2). Note that this value applies only 
to lines forming in the ionized gas. 

The variation in intensity along the amalgamation of 
slit positions 1 and 2 of a variety of interesting line tracers 
is shown in Figs. ^| to [?]. Figure |8| shows intensity varia- 
tions in a number of lines from the IRSPEC data. Figure 
^| plots the intensity profile of selected lines at slit posi- 
tion 3. ^,From Figs. |||7|, we can see that region A coincides 
roughly with a peak in the lines of HI, which have a second 
peak in B, while the emission of molecular hydrogen has a 
strong peak in C and a weaker one in B. Lines of Fell and 
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Fig. 2. J band spectrum obtained towards the three slit sections A, B and C and in the exposition shown in Fig.l. The upper 
spectrum in each panel has been multiplied by 15 to emphasize weak features. 



01 have a sharp peak in B and a secondary one in j4, but 
are absent in region C. In the following, we will use the 
definition A, B and C (see caption to Fig. |^) both to iden- 
tify the peaks in the intensity profiles and with reference 
to Tabic 1. 

We now summarize our results starting with the IR- 
SPEC data which have the advantage that the slit is ori- 
ented perpendicularly to the bar. We then discuss in turn 
the hydrogen and helium recombination lines (which we 
presume form in the ionized gas), the two oxygen transi- 
tions, the collisionally excited iron lines which may form 
close to the ionization front, and the molecular hydrogen 
lines which are thought to form in hot neutral gas close to 
the ionization front. 



3.1. IRS PEC cut 

The observations made using IRSPEC provide a useful 
introduction to the TIRGO results which have a wider 
spectral coverage. Figure || compares profiles in the Brl2 
line from the ionized gas, in the Fell 1.644/xm (4sa 4 D 7/2 - 
3d 7 a 4 F 9 / 2 ) line which traces gas close to the ionization 
front (see below), and in the molecular hydrogen v=l-0 
S(l) line from hot (T> 1000K) molecular gas. Figure || 
demonstrates the fact that the molecular hydrogen peak is 
offset ~16"(0.035 pc) from the ionization front as marked 
by Fell (or by the fall-off in Brl2). The simplest expla- 
nation of this (see Tielens et al. 1993, van der Werf et al. 
1996) is that one is observing an edge-on PDR and that 
the observed offset corresponds to the difference between 
the ionization front where Lyman continuum photons are 
absorbed and the H2 dissociation front where photons ca- 
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Fig. 3. H band spectrum. 



pable of dissociating molecular hydrogen are absorbed. We 
note also that an offset of 16" in the IRSPEC data corre- 
sponds (given the orientation NE-SW of the bar) to 23" 
in the TIRGO slit oriented N-S. 



3.2. HI lines 

In the LONGSP data, we detect many recombination lines 
of atomic hydrogen, 13 lines of the Brackett series (Br7 in 
the K band and 12 lines, from (10-4) to (21-4) in the H 
band), and Pa/3 and Pa7 in the J band. 

We use the 13 lines in the Brackett series to check the 
accuracy of our data. Figure [To] plots the ratio (n-4)/(13- 
4) as a function of the quantum number n and compares 



our results with the prediction of recombination theory 
(Storey & Hummer 1995). The agreement is quite good, 
well within our estimate of 10-20% for the observational 
errors. 

The ratio of Pa/3/Br7 provides a value of the extinc- 
tion Ay~2 mag, assuming the reddening curve of Cardelli 
et al. (1989) and i?=5.5. We detect a slight variation along 
slit 1+2, from 2.3 mag in A to 1.4 mag in C. This vari- 
ation, however, is within our estimated 25% error in this 
line ratio. The value in the CS position is Ay = 1.6 mag. 
Figure [l(] shows as a dashed line the theoretical line ratios 
corrected for Ay— 2 mag. 
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Fig. 4. K band spectrum. 



5.5. #e lines 

One of the aims of our observations was to examine the ex- 
tent to which helium is neutral within the zone of ionized 
gas. Estimates of the helium abundance based upon mea- 
surements of either radio or optical recombination lines 
often assume that the helium and hydrogen Stromgrcn 
spheres are coincident with one another (see e.g. Mezger 
1980). Our profiles along the slit allow us to make a direct 
comparison of the Hel and HI line intensities which can 
then in principle be transformed into the abundance ra- 
tio [ife + ]/[.ff + ] in the immediate vicinity of the ionization 
front of the Bar. 

The chief obstacle in doing this is the uncertainty in 
helium line intensities which results from collisional ex- 
citations from the metastable 2 3 S and 2 : S states. Smits 
(1996) has computed helium line intensities in an approx- 
imation where collisions (and self-absorption) out of the 



metastable levels into n=3 and 4 are neglected although 
collisions between the n=2 levels are considered. We have 
compared our observed intensities with Smits predictions 
for electron density 10 4 cm~ 3 and temperature 10 4 K. 
We normalise for this purpose to the 1.701/im 4 3 D-3 3 P° 
transition which has the same upper level as the 4471 A 
line often used in optical analyses. It is expected that this 
transition (see Osterbrock et al. 1992) is only affected at 
the 1-2 percent level by the collisional effects mentioned 
above and we neglect such effects in the following. One 
finds then that, relative to 4 3 D-3 3 P°, lines such as 5 3 D- 
3 3 P are in good agreement with the Smits predictions but 
4 1 S-3 1 P° and ^S^P at 2.113 and 1.341/xm respectively 
are factors of roughly 3 stronger. The reason for this may 
be the neglect of the collisional effects and trapping dis- 
cussed above (see e.g. Robbins & Bernat 1973; Pcimbcrt 
& Torres- Peimbert 1977). We in any case have assumed 



8 



A. Marconi et al.: NIR Spectra of the Orion Bar 




Fig. 5. Variation with declination offset (relative to declination 
(1950) =—05° 24' 51") of line intensities measured in the amal- 
gamation of slit positions 1 and 2. The lines are H (7-4) (Br7; 
solid) , H 13-4 (dotted), and the He line at 1.701 (im (dashed). 
The vertical scale is arbitrary. We note that regions A, B, and 
C discussed in the text are defined as follows : A, A8 > 16", 
B, -10" < A<5 < 16", C, AS < -10". 

that the 1.701 /im line behaves essentially as predicted by 
the Smits models and can be used to estimate the He + 
abundance. It is natural to compare the He 1.701/^m in- 
tensity with that of the adjacent BrlO line. With the above 
assumptions, we find that : 

F(He, 1.7n)/F(H, 10 - 4) = 3.61 (1) 

In Fig. [ll], we show the profiles along the slit of the 
abundance ratio [i?e + ]/[i? + ] deduced from Eq.(l) as well 
as the ratio of the 2.06 (2 1 P°-2 1 S) to the 1.7/mi lines. The 
Br7 profile along the slit is shown for comparison. We 
derive an abundance ratio of 0.093 ± 0.005 over region A 
consistent with He abundance estimates from other au- 
thors (see e.g. Baldwin et al.1991 who find 0.088 based 
on optical measurements and the review of Mezger 1980 
(his Fig. 5) who shows that radio estimates at positions 
less than 100" from 6 1 C Ori are in the range 0.083-0.09). 
It appears that at the position of our slit, the helium and 
hydrogen Stromgren spheres are close to being coincident. 
We note with interest the decrease of [i?e + ]/[iJ + ] to val- 
ues of ~ 0.075 in the southern part of the slit or effectively 
in zone C, where we seem to observe different ionization 



Fig. 6. Variation with declination offset of H (7-4) (dashed 
line) and H2(1-0)S(1) (solid line), measured in the amalgama- 
tion of slit positions 1 and 2. The regions defined as A, B, and 
C are shown. 

conditions in this (presumably) foreground ionized mate- 
rial. One notes also the fact that the degree of "enhance- 
ment" of the 2.06 micron line seems to increase slightly 
(by a factor of 1.3) in zone C. 

3.4- OI lines 

We detect two OI lines in the J band, the 2p 3 3D-2p 3 3P at 
1.129/zm and the 2p 3 4S-2p 3 3P at 1.317/^m, with compara- 
ble intensity. These lines are produced in the neutral gas 
by excitation of OI to the upper level of the transitions by 
UV photons, at 1027 and 1040 A respectively, followed by 
radiative decay. The fact that the two lines have similar 
intensity (though not at the CS peak, see Fig. 2) suggests 
that the contribution of Ly— (3 photons to the excitation 
of the 2p 3 3D level (the upper level of the 1027 A line) is 
negligible. That the excitation mechanism is fluorescence 
is confirmed by the spatial variation of the line intensities, 
which peak at the edge of the ionized region as marked by 
the HI recombination line emission (cf. Fig. 0). In other 
words, the oxygen lines are an excellent marker of the ion- 
ization front. More precisely, one can say that due to the 
rapid charge exchange of OI with H + and the fact that 
the hydrogen and oxygen ionization potentials are close 
to being identical, the oxygen lines trace neutral gas close 
to the ionization front. 
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Fig. 7. Top panel: variation with declination offset of H (7-4) 
(dotted line), H2 f-OS(l) (dashed line) and Fell f.644 ^m(solid 
line). Bottom panel: variation with declination offset of H (7-4) 
(dotted line), H2 f-OS(l) (dashed line) and OI 1.317 ^m(solid 
line). In both cases, the variations are measured in the amal- 
gamation of slit positions 1 and 2. 



It is possible to use the 01 lines to derive a measure of 
the UV radiation field Go at the edge of the bar. If fluores- 
cence is the dominant excitation mechanism, the number 
of photons emitted in the IR line equals the number of 
photons absorbed by the UV line. Since the optical depth 
of the UV lines is always very large (to ~ 5 V/10 18 for the 
1040 A line and - 10 V/10 18 fo the 1027 A where N is the 
atomic H column density and we have assumed a velocity 
dispersion Av=3 km s" 1 and 0/11=6 x 10~ 4 ), the number 
of absorbed UV photons is proportional to the line equiv- 
alent width in the "flat" portion of the curve of growth. 
This is given by : W x /X uv ~ 1.2Av/cF{t ) ~ 3.6 x 10~ 5 
(Spitzcr 1978, p53, where To is the UV line center optical 
depth and F(to) ~ 3 for optical depths of order 1000). 
Since both UV lines are in fact triplets with separation 
larger than W\ , the total number of UV photons absorbed 
and re-emitted in each IR line is given by: 



I(IR) (erg cm s Hz ) (2) 



ju v 



cW x 



where I(IR) is the observed intensity of the IR line in 
erg cm~ 2 s _1 sr _1 and t is the angle between the PDR 
and the line of sight (0 t = 90deg in a face-on PDR; see 



Fig. 8. IRSPEC cuts: HI (12-4) (dashed line), H2 1-0S(1) (dot- 
ted line) and Fell 1.644 /im (solid line). 



Appendix). Note that the three components of the IR lines 
are not resolved in our spectra. 

Equation 2 gives I^ v ~ 1.2 x 10- 13 sin(9 t for 1(1 R) = 
2.6 x 10~ 4 erg cm~ 2 s _1 sr _1 , as observed in the main 
peak of the 1.317 /im line. Here, we have corrected for 
2 magnitudes of visual extinction (see Sect. 3.2). The in- 
ferred UV intensity can be compared to the flux from 6 1 C 
Ori at the projected distance of the bar (I p ~4x 10~ 14 
erg cm-2 s" 1 Hz" 1 for T*=40000K, L*=2.5xl0 5 L ). 
The main uncertainty is the appropriate value for t but 
the Orion Bar is known to be close to edge on (see e.g 
Hogcrhcijdc ct al. 1995). A plausible value is sin# t <~ 0.2 
(0 t ~ 10 — 15deg; see Sect. 3.6), and one finds then that 
the physical distance of the bar from 1 C Ori is about 
(I p /I^ v ) ' 5 or 1.3 times the projected distance. 

The UV intensity can be expressed relative to the inter- 
stellar diffuse field taken here to be 3 x 10 photons cm~ 2 
s _1 . The normalized UV intensity is then 1.3 x 10 5 sin0 t 
or Go ~ 2.6 x 10 4 (for sin# t =0.2), similar to the value 
used by Tielens et al. (1993). 

The second, weaker peak of emission A at AS ^23 
arcsec has an OI intensity about two times lower than the 
main peak. This may be due to a different orientation of 
the front with respect to the line of sight. The OI intensity 
on the CS peak is about 1 /4 that of the main peak which 
again may be due to an orientation effect although it could 
also imply dust extinction between the CS peak and the 
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Fig. 9. Variation with declination offset along the slit cen- 
tered on the CS peak of selected lines: HI (7-4) (dashed line, 
top panel), H2 1-0S(1) (solid line, top panel and dotted line, 
bottom panel), and Fell 1.257 ^m(solid line, bottom panel). 



Trapezium. In general, our 01 results imply Go values at 
the ionization front in the range 6000 — 3 x 10 4 . 

3.5. Iron lines 

Model calculations suggest (see Baldwin et al. 1991, Rubin 
et al. 1991, Osterbrock et al. 1992) that iron is mainly in 
the form FelV in the Orion nebula. It follows that one ex- 
pects to see Fell and Felll emission predominantly at the 
edge of the Stromgren sphere close to the ionization front. 
In Fig. 0, we show profiles along our amalgamated slit of 
the 1.644 (im [Fell] transition compared with the 1.317/im 
OI line discussed in the previous section. One sees that the 
two lines show rather similar behavior with a peak slightly 
offset from one of the maxima observed in the H recombi- 
nation lines. It seems plausible that this approximate co- 
incidence denotes the presence of an ionization front and 
the OI data which we discussed above confirm this idea. 
It is notable also that the iron lines show no evidence for a 
coincidence with the molecular hydrogen 1-0 S(l) Cpeak 
at A<5= -27" which is also shown on Fig. [?]. In fact, our 
data suggest that both the Fell and Felll emission lines 
form in ionized (or partially ionized but not neutral) gas 
close to the ionization front. It is worth noting here that 
extinction estimates which we have made using the ratio 
of the 1.644 4sa 4 D 7/2 -3<i 7 a 4 F 9/ 2 to the 1.26^m 4sa 4 D 7/2 - 



Fig. 10. Ratio of the H-band hydrogen recombination lines 
(n-4)/(13-4) as a function of the quantum number n. Filled 
circles refer to position A, triangles to B, squares to C, and 
diamonds to position CS (the CS peak). The line (18-4) is 
blended with an [Fell] line. The solid line shows the predictions 
of recombination theory (Case B). The dashed line shows the 
Case B ratios corrected for a reddening of Av~2 mag. 



4sa 6 Dg/ 2 transitions is consistent with a visual extinction 
of 2.7 ± 0.9 magnitudes at all positions. 

We can estimate the conditions required to explain the 
relative intensities of the Fell lines. Fell line ratios can 
be used as indicators of electron density (see Oliva et al. 
1990, Pradhan & Zhang 1993). Our most useful indicator 
appears to be the ratio of the 1.600 4sa 4 D 3 / 2 -3(i 7 a 4 F7/ 2 
line intensity to that of the 1.644 4sa 4 D7/ 2 -3<i 7 a 4 Fg/2- 
We find that this ratio varies in the range 0.06-0.1 (see 
table 1) over the regions covered by our slit. Based on 
the collisional rates of Pradhan & Zhang(1993) (see also 
Oliva et al. 1990), we conclude that this corresponds to an 
electron density n e of 4000 cm~ 3 in the region ^4), 6000 
cm~ 3 in region C, and 3000 cm -3 in region B. At the CS 
peak position, the observed 1.600/1.644 ratio is 0.12 cor- 
responding to n e = 10 4 cm~ 3 . Thus, we can exclude high 
density clumps with n e = 10 6 cm" 3 of the type discussed 
by Bautista et al.(1994). On the contrary, the Fell data 
seem consistent with ionized gas of electron density ~ 10 4 
cm -3 or less in the vicinity of the ionization front. It is 
worth stressing here that this line ratio converges to the 
LTE value at electron densities of roughly 10 5 cm~ 3 and 
so the lack of high density clumps is significant. More- 
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Fig. 11. The figure shows the He abundance [He + ]/[H + ] and 
the ratio of the two Hel lines at 2.058 and 1.701 /im as a 
function of the declination along the amalgamation of slits 1+2, 
compared with the Br7 profile. The data have been smoothed 
over three pixels. 



over, the ionization degree in the layer where the Fell lines 
arc formed is unlikely to be much smaller than unity and 
hence the hydrogen density is also likely to be of order 10 4 
cm" 3 . 

3.6. H2 lines 

We measured the intensity of 8 H2 lines in the K band, cov- 
ering an excitation range from 6472 to 18089 K (see Table 
1). Using the data from table 1 and transition probabilities 
from Turner et al. (1977), we have determined upper level 
column densities at positions A, B, C and CS. In Fig. |l2|, 
we plot the column densities per sub-level derived in this 
manner against excitation energy. One sees that there are 
clear departures from an LTE distribution suggesting ei- 
ther that fluorescence is playing a role in determining level 
populations or that there is a sharp gradient of tempera- 
ture along the lines of sight sampled. The "best fit tem- 
peratures" derived from fitting a Boltzmann population 
distribution to the data in Fig. 12 are moreover rather 
high with values ranging from ^2500 K in region C to 
3000 K in region A. 



The high excitation temperatures as well as the (prob- 
able) detection of lines from levels as high as v=3 suggest 
that we are detecting extended fluorescent emission (seen 
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Fig. 12. Column density per sub-level against excitation en- 
ergy of the level at positions A, B, Cand on the CS Peak. The 
dotted lines show the best fit to a single temperature popula- 
tion. The two v=3 transitions are doubtful and have not been 
considered on the fit. Filled squares represent ortho transitions 
and circles para transitions. 



on larger scales by Luhman & Jaffe, 1996) in addition to 
a "thermal" layer. 

The evidence that some fluorescent emission is present 
is strengthened by the behaviour of the intensity ratio 
of the v=l-0 S(lj (2.12/xm) and v=2-l S(l) (2.25 /jm) 
H2 lines along our amalgamated slit (A,B,C), shown in 
Fig. |l3|. In a pure fluorescent model, this ratio is predicted 
to be ~2, whereas an admixture of collisional excitation 
leads to higher values (12 for pure thermal emission at 
2000 K). Our results are consistent with those of van der 
Wcrf et al. (1996) in that at the main molecular hydrogen 
intensity peak C, the ratio is ~ 6-7. A similar value is also 
found at the secondary peak B, which is about coincident 
with the position of the ionization front as traced by the 
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Fig. 13. Ratio of the 1-0S(1) to 2-lS(l) H 2 lines as a func- 
tion of the declination offset along our amalgamation of slits 
1+2 (solid line). To compute the ratio, the profiles have been 
smoothed over three pixels. The horizontal bar shows the value 
averaged over declination offset 0-50 arcsec. We show for com- 
parison the profile of the 1-0S(1) line (dotted line). 

01 lines. Over the rest of the slit the ratio is <4. Thus the 
fraction of fluorescent emission is smaller at the peaks of 
v=l-0 S(l) emission. 

These measurements, as well as the observed intensity 
of the v=l-0 S(l) line, can be compared to the PDR model 
calculations of Hollenbach & Natta (1995) (at steady 
state) . Figure [l5| shows the predicted intensity of the 1-0 
S(l) line in a face-on PDR as a function of the density 
for different values of Go, and in Panel (2) the ratio of 
the 1-0 S(l) to the 2-1 S(l) line. In the main molecular 
peak C, we observe a peak intensity of the 1-0 S(l) line 
of ~ 2.3 x 10 -4 erg cm~ 2 s^ 1 sr -1 , and a ratio of ~7. If 
we assume that the H2 lines are produced in a PDR with 
Go ~ 3 x 10 4 (see Sect. 3.3), this corresponds to a density 
6 x 10 4 cm" 3 . 

These results assume a face-on PDR and there is good 
reason to sup pose th at the Orion Bar is seen edge on (see 
Jansen et al. |l995b| for a discussion of the geometry). The 
effect of a slant of the Bar on the H2 lines is complicated 
by the effects of dust extinction and a brief discussion is 
given in the appendix. The results are different for the vi- 
brationally excited H2 lines at 2/j.m and for lines at longer 
wavelengths, for which the extinction is negligible. Par- 



Fig. 14. Same as Fig. 13 for the CS position. 

mar et al. (1991) have observed the J=3-l (17/im) and 
4-2 (12.3/im) v=0 H2 lines toward the bar and find that 
their observations are compatible with a hydrogen column 
density of 3 x 10 21 cm -2 at a temperature of order 500 
K. Figure |l^ shows in Panel (3) the model-predicted ra- 
tio for the two lines observed by Parmar et al. and in 
Panel (4) the ratio of the v=2-l S(l) to the J=3-l (17^m) 
lines. These ratios can be well reproduced by a model with 
G ~ 3 x 10 4 and n ~ 6 x 10 4 cm having a moderate 
enhancement of the intensity of the 12 and 17 jum lines 
due to an inclination of the Bar with respect to the line 
of sight 6 t ~10deg. These values are also consistent with 
the inclination required to explain the intensity of the OI 
lines, discussed in Sect. 3.4. 

The density in the other regions where molecular lines 
are measured can be estimated in a similar way from Fig. 

In the secondary molecular peak B, assuming the same 
inclination t ~ lOdeg, we obtain n ~ 4 x 10 4 cm~ 3 . 
These densities are a factor of ~5-10 larger than that of 
the ionized gas (see Sect. 3.5) and roughly consistent with 
the density required to explain the stratification of the bar 
(i.e the offsets between cold molecular gas and ionization 
front, see Tielens et al. 1993), which is of order 5 x 10 4 
cm~ 3 . 

Similar densities can be derived from our data at slit 
position 3 , the "CS peak" {CS). We show the 1-0 S(l) in- 
tensity in Fig. ^| and the v= 1-0/2-1 line ratio in Fig. [l4| 
We see a peak of H2 emission at Ad ~ 49"with an intensity 
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Fig. 15. Panel (1): model predictions for the intensity of the H2 1-0 S(l) line as a function of hydrogen number density for 
different values of Go- Panel (2): ratio of the v=l-0 S(l) line at 2.12 /im to the v=2-l S(l) line at 2.24 /im as a function of 
density. Panel (3): ratio of the v=0,J=4-2 line at 12 /im to the v=0,J=3-l at 17 /im. Panel (4): ratio of the v=l-0 S(l) line to 
the v=0,J=4-2. Each curve is labelled by the logarithm of Go- Triangles indicate Go = 10 5 , squares with Go = 10 4 ' 5 , circles 
with Go = 10 4 , and crosses with Go = 10 3 . The observed value at peak C is shown by the arrow in each panel. 



of ~ 1.6 x 10~ 4 erg cm -2 s" 1 sr -1 , and a broader emission 
between 0"and 30", with a peak intensity of ~ 8 x 1CP 5 erg 
cm~ 2 s _1 sr _1 . In both cases the ratio of the two H2 lines 
is about 5 implying densities <10 5 cm~ 3 for Go ~ 10 4 . 
Neither of the peaks coincide with the "ionization front" 
as defined by Fell (Fig. ||). One would expect such a coin- 
cidence if the densities were considerably above 10 5 cm -3 
as implied by the CS data (see van der Werf et al.1996). 
Thus the general conclusion is that the PDR as seen in 
molecular hydrogen appears to be at densities below 10 5 
cm -3 in contrast to the millimeter results which suggest 
the existence of clumps at densities well above 10 5 cm~ 3 . 

A caveat to much of the above discussion is that com- 
parisons which we have made between the predictions of 
the Hollcnbach & Natta code used by us and other results 
(in particular, the results of Storzer & Hollcnbach 1997) 
shows that there can be substantial differences in both the 



H2 line intensities and in some line ratios. The main rea- 
son for this is the extreme sensitivity of H2 line intensities 
to the temperature structure although there are more mi- 
nor effects which also play a role (note for example that 
the Hollcnbach Natta models neglect the effects of the car- 
bon oxygen chemistry upon H2). This has the consequence 
that rather minor errors in the treatment of thermal equi- 
librium can affect our conclusions. In this respect, mid 
infrared intensity ratios such as the 12/17 micron ratio 
measured by Parmar et al. are important because they 
afford a direct measure of temperature. Thus if the tem- 
perature structure can be adjusted to fit the mid infrared 
observations, one may have some confidence in the predic- 
tions for other lines. For the moment, we conclude that the 
fact that we have been able to fit the 12/17 ratio suggests 
that the model which wc have used is reliable. 
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4. Discussion 

Our analysis of molecular hydrogen in the previous sec- 
tion has assumed implicitly that the H2 emission comes 
at all positions from material heated by the stellar radia- 
tion field (PDR) . The two peaks C and B we observe must 
then come from two different structures. In fact, there are 
indications in the molecular hydrogen images of van der 
Werf et al. that in H 2 lines there are two bars. One of these 
(which one might call the main H2 bar) corresponds to the 
feature seen in our Fig. || at position C (AS = —24"). The 
"second H2 bar" (less well defined in the van der Werf et al. 
images) is close in projection to the main ionization front 
and corresponds to the peak in Fig. || at AS=-2" (roughly 
coincident with the 01 peak). We propose that these two 
bars are separated along the line of sight and that the shift 
in position (15" accounting for our slit orientation) is due 
to a tilt in the bar of ~10 degrees as discussed above. 
Thus, the bar is split along its length (essentially along 
the line of sight) into two sections separated by 0.2-0.3 
parsec. Each half of the bar in this scenario has a length 
of order 0.1 parsec and thus the total length along the 
line of sight is of order 0.3-0.4 parsec or somewhat smaller 
than proposed by Jansen et al. ( 1995b| ). 

There are problems with this model however. The prin- 
cipal difficulty is that our data do not show evidence for 
ionization front indicators (OI and Fell) roughly coinci- 
dent (one expects a shift of order 6 arc second for a den- 
sity of 5 10 4 cm -3 ) with the main (AS = —24" on our slit) 
molecular hydrogen peak. One explanation for this might 
be that the density in the layer of gas between the ion- 
ization and H2 dissociation front is a factor of 2-3 lower 
(of order 2 x 10 4 cm -3 ) and thus that the main H2 peak 
is shifted 15" (20"in our NS oriented slit) relative to the 
main ionization front (i.e., peak C in H2 corresponds to 
peak B in OI ). Equally, the ionization front emission seen 
in Fig. [?] at A<5 =23" might correspond to the H2 emis- 
sion at AS = 3". This lower density between ionization 
and photo-dissociation fronts implies a density increase 
between the atomic and molecular regions by at least a 
factor of 4 since the general bar stratification requires an 
average density of at least 5 x 10 4 cm -3 in order to account 
for the observed offsets of ionization front and molecular 
lines (see also Wyrowski et al. 1997). 

Thus, one possible interpretation (see also Simon et 
al. 1997| ) of the H 2 data is that there is a sharp density 
gradient perpendicular to the line of sight such that the 
molecular gas has much higher density than the partially 
ionized atomic medium adjacent to it. This has the at- 
tractive feature that it helps explain one of the puzzles 
concerning the bar which is the discrepancy (more than 
an order of magnitude) betwee n the hydrogen column den- 
sity derived by Parmar et al. ( 1991 ) and that inferred by 
Hogerheijde et al. (|l995|) from their C 18 data. The Par- 
mar et al. data refer essentially to the main H2 bar whereas 
Hogerheijde et al. preferentially sample the fully molecu- 



lar gas to the SE. A density gradient may also cause the 
offset between the ionization front and the main molecular 
hydrogen peak to be larger than that estimated using a 
homogeneous model whereas molecular hydrogen and car- 
bon radio recombination lines would become closer to one 
another. 

This last aspect is particularly interesting in view of 
the coincidence found by Wyrowski et al. (1997) between 
the bar seen in C91a emission and the main H2 bar. The 
difficulty in explaining this result stems from the fact that 
one expects the H2 emission to come from gas with tem- 
perature above 2000 K while the carbon line is thought to 
be formed at temperatures which are considerably lower. 
In fact, there is a firm upper limit of 1600 K on the temper- 
ature of the gas responsible for the C91a emission (based 
on the line width). The proposed density gradient dis- 
cussed above may cause the offset between the photodis- 
sociation front (i.e H2 emission) and carbon line emission 
to diminish. 

A different interpretation of our observations is in prin- 
ciple possible. Our main molecular hydrogen peak C at 
AS = —24" could be produced by a low velocity shock pre- 
ceding the ionization front. This would explain the lack of 
ionization front indicators coincident with molecular hy- 
drogen. However, the kinematics of the emission seen in 
C91a by Wyrowski et al. are difficult to explain in this 
scenario . One needs a shock velocity of at least 3 kms -1 
to excite H2 (see Tielens et al. 1993) and then the ob- 
served C91a line widths become difficult to understand. 

It is intriguing that our H2 observations do not show 
any evidence of high density gas. This is in contrast with 
the fact that the molecular line data (e.g. Tauber et al. 
1995, Simon et al. |1997|, van der Werf et al. 1996) give 



evidence for a considerable fraction of the gas being in 
clumps with density 3> 10 5 cm -3 . Such clumps can be 
expected to affect our results because most of the lines 
observed by us are sensitive to high emission measure and 
high density. High density PDRs are expected to be hot- 
ter and hence considerably more intense in H 2 v=2-l and 
1-0 emission than lower-density PDRs (cf. Fig. [IB]). Nev- 
ertheless, our data show no evidence for gas with density 
above 10 5 cm -3 even towards regions where Simon et al. 
(see also van der Werf et al.) estimate molecular hydro- 
gen densities of order 2 10 5 cm -3 . We see no reason on 
the other hand why clumps should dissipate on a short 
timescale when traversing the H2 dissociation front and 
suggest therefore that clumps, while possibly present, are 
a secondary phenomenon. 

More important in our opinion is the density gradient 
mentioned above. It is worth noting that in the scenario 
which we are advocating, the thermal pressure may be 
constant (at a value of the order of 10 8 cm _3 K) along a 
line perpendicular to the bar and we conclude that isobaric 
models of the Bar are worth examining. This incidentally 
would suggest relatively low values for the magnetic pres- 
sure and hence magnetic field (below 0.5 mG). 
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5. Conclusions 



This study has presented NIR slit spectra of the Orion Bar 
region. Our main result is based on the molecular hydro- 
gen line intensities and is that the densities derived from 
these tracers are of order 3 — 6 x 10 4 cm -3 and thus consis- 
tent with estimates of the mean density derived fro m the 
observed stratification of the bar (e.g. Tielens et al. 1993| 
Wyrowski et al. 1997). Comparison with the longer wave- 
length H2 data of Parmar et al.(1991) implies a tilt for 
the bar relative to the line of sight of ~10 degrees. Our 
data suggest also that the bar may be split into two por- 
tions along the line of sight which are separated by 0.2- 
0.3 parsec. It seems plausible that the density is lower in 
the "atomic" region (perhaps 2 x 10 4 cm~ 3 ) than in the 
molecular gas (of order 10 5 cm -3 ). We conclude that mod- 
els with constant thermal pressure should be examined in 
future studies of the bar region. 

We have also derived densities for the ionized gas in 
the vicinity of the ionization front using [Fell] line ratios 
and find values of order 10 4 cm -3 . This together with the 
molecular hydrogen data has convinced us that high den- 
sity neutral clumps play a rather minor role in determining 
the observed characteristics of the bar. 

A by-product of these observations was that we were 
able to use the observed 01 line at 1.317^im as an esti- 
mator for the ultraviolet radiation field incident upon the 
bar. We estimate the normalised FUV intensity Go on the 
bar to be 0.6 — 3.0 x 10 4 using this tracer. 

Finally, we have used the 1.701/im line of He to ex- 
amine the degree of coincidence of helium and hydrogen 
Stromgren spheres. To within our errors, we find that He + 
and H + coexist and hence that He abundance estimates 
using these tracers should be reliable. 
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A. The effects of inclination on the H 2 lines 

The intensity of a given line as a function of the angle 9t 
can be approximately estimated from the expression (cf. 
Fig. O): 
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Fig. 16. Sketch of the geometry for molecular hydrogen emis- 
sion discussed in the Appendix. The symbols used in Eq.(Al) 
are defined. The observer sees the PDR at an angle 9t- The 
layer of H2 emission has continuum optical depth to and the 
foreground atomic layer optical depth ro,„. 



frequency across the layer of the PDR which contributes 
to the H2 emission, r„ the optical depth of the HI re- 
gion along the line of sight, tq and To, ra are the optical 
depths in the H2 and HI regions, respectively, in the di- 
rection perpendicular to the PDR. To first approximation, 
t ~ To/sin# t , r n ~ To, rl /sin# t . For 9 t <l/L, the line of 
sight does not intercept the neutral region of the PDR 
(t„ ^0) and the line intensity tends to I/I ~ 1/t . When 
to = To, n = 0, then I/I = l/sin6l t . 

Figure [l?] shows the ratio I/Iq for the 1-0 S(l) line 
(solid curve), which has been derived using tq ~ 1 x 
A\/A v , t ,„ ~ 1 x A\/A v , A\/A V =0A4:5. Here we assume, 
based upon models, 1 magnitude of visual extinction in 
the H2 emitting layer and 1 magnitude of extinction in 
the foreground HI layer. For the 1-0 S(l) line, I/Iq ~ 1 
with an accuracy better than 30% for 8 t > 8deg . The 
figure shows also the variation with 9 t of the ratio of the 
1-0 S(l) line at 2.12 /j,m to the v=0,J=3-l 17^m line (for 
which To ~ To )Tl ~ 0), with respect to the same ratio for a 
face-on PDR. 



/o/r (l - e~ T )e 



(Al) 



where Iq is the face-on intensity (9 t = 90) ignoring ex- 
tinction due to dust, r the dust optical depth at the line 
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Fig. 17. The solid line shows, as a function of the inclination 
angle 9t, the intensity of the 1-0 S(l) line normalized to the 
value in a face-on PDR where dust extinction is neglected. The 
dashed line shows the ratio of the two lines 1-0 S(l) and v=0, 
J=3-l at 17 /jm, normalized in an analogous faction. 
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